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I. INTRODUCTION 



In order to explain the current accelerated expansion of the Universe, we need to modify 
the Einstein equation. There are two ways to modify the equation. One way is the modi- 
fication of the right hand side of the Einstein equation by including unknown cosmological 



fluid, called "dar 
/c-essence model 



t energy" . A typical example of the dark energy is the ACDM model and 
l sj] is also one of the dark energy models. The /c-essence model had been 



known as the model which causes the accelerated expansion of the Universe only by the 
kinetic terms of a scalar field. However, the /c-essence model is currently recognized as a 
general single scalar field model which includes quintessence model {4- 7], ghost condensate 
model [8j, tachyon model UQ and so on. Another way to modify the Einstein equation 



is the modification of the left hand side of the Einstein equation, which is callec 
gravity". For example, scalar-tensor theory 
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151 ] and F(R) gravity model 



"modified 
are 
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known as the modified gravity models. Only from the data about the expansion rate of the 
Universe which is given by the observations of supernovae of type la, we cannot determine 
which model could be real one among such many models. One way to solve the problem is 
to consider the evolution of the matter density perturbation. The behavior of the evolution 
of the matter density perturbation is usually obtained by solving the equations of the linear 
perturbations, consisting of the metric perturbation and the perturbation of the energy mo- 
mentum tensor. It is difficult, however, to solve the equations because scalar field models 
have more degrees of freedom than those in the ACDM model. Hence, the remaining way 
to evaluate the evolution of matter density perturbation has been using computers to solve 
the equations numerically under special conditions or to apply subhorizon approximation to 
the equations. 

In this paper, we will see that by linearizing equations without any other approximation, 
we find the behavior of the evolution of the matter density perturbation in Sec. IIHI We 
will also see that the matter density perturbation shows the oscillation whose frequency 
is proportional to the sound speed in /c-essence model at the leading order of small scale 
and the effective growth factor of it is determined up to by the second derivatives respect 
to X = —d fJ- (j)d fl (j)/2 of the Lagrangian density K((p,X) and and their time derivatives. 
A clear observational difference between the ACDM and /c-essence model could be also 
developed by the sound speed in /c-essence model, which is also shown in Sec. IIIII The last 
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section is devoted to the conclusions and discussions. We use units of &b = c = H = 1 
and denote the gravitational constant 8irG by k 2 = 8n/Mp 2 with the Planck mass of 
Mp, = G~ 1/2 = 1.2 x 10 19 GeV. 



II. COSMOLOGICAL PERTURBATIONS 

A. Equations of linear perturbations 

We start with the following action of /c-essence model: 

S = J d i x^^-K( ( p,X) + L mattCT ^, X = ~d^d^. (2.1) 

Here is a scalar field and L mattCT expresses the Lagrangian density of matters. Following 
the variational principle, by the functional differentiation with respect to metric g^ix), we 
obtain the Einstein equations, 

-i^ = - K(4>, X)g^ - K x d^d u cf> + Tj£> , (2.2) 

where K x expresses the partial derivative of K((f),X) with respect to X and is the 
energy-momentum tensor of the matters. On the other hand, the equation of motion for the 
scalar field <j) is given by 

K 4 , + V (t (K tX & i <l>) = 0. (2.3) 

When the Universe is described by the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
space-time with the flat spacial part, whose metric is given by ds 2 = —dt 2 + 
J2i j-i 2 3 5 i ja 2 {t)dx' l dx^ , Eqs. (12 .2p and (I2.3P have the following forms: 
2>H 2 

=K - 4> 2 K, X + p matter , (2.4) 

K 

- - <j) 2 K tX + (1 + w)p 

matter • 

(2.5) 



K 2 



=3H(t)K x + <P{K )X + <j> 2 K xx ) - K+ + <f) 2 K^ , (2.6) 

where w = p m attcr/Pmatter is the equation of state (EoS) parameter and H{t) = a(t)/a(t). 
The energy momentum tensor is treated as that of perfect fluid, here. Then, the continuity 
equation for the matters is given by 

Pmatter + 3(1 + w)Hp maitteT = . (2.7) 



In the following, we abbreviate the suffix "matter" for the quantities corresponding to matter, 
e.g., we write energy density p ma tter by p, etc. We now consider the perturbed scalar field 
and the metric around the background, — > + 5<f>, X — > X + 5X, and g^ u — > g^ v + Sg^, 
to analyze the time evolution of the matter density perturbation. Then Eqs. (12. 2p and (12. 3p 
give the following equations of linear perturbations, 

- ^6K >x d^<j> - l -K, x (d^d»<P + d^d v 8cj>) + 5T™\ (2.8) 

=K M 5<t> + K^ X 5X + V^6K x d»<f> + K, x d a <J)5g a » + K^Stfy 

+ 5Y^K^. (2.9) 

Here, we use the metric of FLRW space-time with the Newtonian gauge, 

ds 2 = (-l + 2$)rft 2 + S ij a 2 (t)(l + 2^)dx i dx j . (2.10) 

ij=l,2,3 

Then, Eq. (I2.8P has the following forms: 

k 2 ( 
-6H 2 $ - 2—^ - 6H^> =k 2 { i-K x + <p 2 K xx )U 2 <& + 050) 
er I 

- {K+ - <P 2 K X ^ - 20A>W + 2<S><p 2 K x - 8p} , 

(2.11) 

2di(H$ + if) =k 2 {j>K x di5<f> + (p + p)S Ul } , (2.12) 
a- 2 ^($-M/)=0, (i^j), (2.13) 

nf_QA_ 2Hdo _ Ak _ 6H2 \ ^_f^ + m + 2dodo + 6Hd \ ^ 

\a z a z J \a z a z J 

=k 2 {-if x (0 2 $ + 050) - K^Scj) + 5p} , 

(not summed with respect to i), (2-14) 

where Eqs. (I2.1ip -( l2.14p express the (0,0), (0,z), for i ^ j, and for i — j 

components of the Einstein equations, respectively. In Eqs. (12.1 ip and (12.141) . k expresses 
the wave number which appears from the derivative with respect to the spacial coordinates 
(k 2 = —djdj) by the Fourier transformation. 

Since we have treated energy momentum tensor as perfect fluid T^ v = pg^ v + (p + p)u^u u , 
where it, = and u = —1, the perturbation of the energy momentum tensor is given as 
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follows, 

ST° = -Sp, (2.15) 

8T° i ^(p + p)5u i: (2.16) 

6T = -a- 2 (p + p)5 Ui , (2.17) 

S^^Sp, (2.18) 

where we have used the relation 5uq = 5goo/2, which is obtained from the condition 
g^u^Uy = —1. We now decompose 5ui as d{5u + Suf, where 5ui satisfies the transverse 
condition d l 5u( = 0. When we consider the scalar perturbation, we put 5u( = and 
keep only the scalar part 8u. Then, 5ui is expressed as 5ui = di5u. On the other hand, 
perturbative equation of the scalar field ( 12. 9ft has the following form: 

(<j> 2 K >xx + K >x )5& = - 3<j>K iX i> + (-0^,x«9 o - <j> 3 K xx d - QH<j)K,x - 20K, X 

- 2j>K iX - 3H^K >XX - 3<p 2 4>K tX x ~ <p 3 K >XX + 2 ^x)$ 

- (3HK, X + K x + 3H<p 2 K, XX + 2^K XX + <p 2 K X x)S<P 
-K jX ^ - 3H<j>K iX<t> - <j)K^ ~ 4>K,Xcf> + 50 . (2.19) 

The perturbative equation of matter is given by linearizing the continuity equation V^T^ = 
0: 

5p + 3H(5p + 6p) + a- 2 di{(p + p)5 Ul } + 3^(p +p) = , (2.20) 
a~ 3 d {a 3 (p + p)8u i ) + didp - (p + p)di<& = . (2.21) 

In the following, we consider the radiation dominant era, the matter dominant era, and the 
subsequent era and therefore we may assume w = constant. In addition, we assume the 
sound speed c 2 = Sp/Sp behaves as constant. 



B. Subhorizon approximation 

By using Eq. ( 12 . 1 3j) . we find we can choose \P = which is assumed in the following. 
Subhorizon approximation is the approximation that H 2 , HOq, d 2 <C k 2 /a 2 . If we apply 
the subhorizon approximation in case of p = 0, c 2 = 5p/5p = i.e. 5p = 0, we obtain the 
following relation from Eq. (12. 19ft . 

> . (2.22) 
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This relation implies we can neglect S(j) in Eq. (12. lip and therefore we obtain the approxi- 
mated equation from Eq. (12.111) : 



The above equation (I2.23|) is identical with the equation of the Newton gravity. We can, 
then, obtain the evolution of matter perturbation by using Eqs. (12.221) . (12.231) . ( 12.201) . and 
(EOT]) . Differentiating Eq. ( I2~20l) with respect to N = In a and substituting Eq. ( |2~2TT) . we 
obtain the following differential equation, 



where 5 = Sp/p and f2 m is the fractional density of matter, defined by K 2 p/(3H 2 ). w c g is 
an effective EoS parameter defined by w c g = — 1 — 2H'/3H. This equation (12.241) does not 
include the contribution from the scalar field <fi explicitly, so that Eq. ( I2.24p is not changed 
from the equation given in the ACDM model. Thus, it might be considered that there is 
no contribution for the matter density perturbation explicitly from the scaler field when 
scaler field is only minimally coupled to the gravitation. This is, however, incorrect. Because 
the hypothesis we are able to use the subhorizon approximation is too strong to obtain a 
correct answer. The subhorizon approximation is decomposed into two approximations, that 
is, small scale approximation k 2 /a 2 ^> H 2 and the assumption that the magnitude of the 
time variation of gravitational potential and the perturbation of the scalar field are almost 
that of the Hubble rate such as do^ ~ Hty, and the latter approximation can be only used 
in the case that the propagating speed of the perturbation mediated by the scalar field is 
negligible compared with the light speed. The behavior of the matter density perturbation 
in case that the contribution from the scalar field is not negligible in the subhorizon scale 
will be investigated in the next section. 

III. PERTURB ATIVE EQUATIONS WITHOUT APPROXIMATIONS 

In this section, we derive the differential equation of the matter perturbation 5, without 
using any approximations except for using the simple Newtonian gauge metric perturbations 
and the perfect fluid approximation. 
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\l> ~ 5p. 



(2.23) 




(2.24) 
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A. Linear differential equations 



First, we rewrite the differential equations in the linearized forms. After that, we will 
obtain a higher order differential equation for the quantity 5 = Sp/p in the next subsection. 
In the following, we again assume \& = $ from Eq. ( I2.13p . Then, Eq. (12. lip yields, 



* - P 5p - Pi* - P 2 5X - P 3 c 
- 2 k 2 







K 



Po= GH> Pl 



-H 



3a 2 H 



,P2 



K 



6H 



K, x + rK iXX 



6H 



K,-6 2 K 



X<f> 



Eqs. (EUD and flXT]) give, 



5X - E 5p - E^ - E 2 6X - E 3 



E n =- 



K 



x 



2H Kx + ^Kjcx 



Ei =4>{K. 



x 



xx , 



M3 H + 



k 7 



3a 2 H 



4>K 



x 



K 



.2 



d 



E 2 =—6 2 K x - ^ln 
2fT * dt 

E 3 =-^(K,x + ^K tX xr 



(K,x + <p 2 K t xx] 



K 

2H 



—<j>K,x(K,<f> - 2 A> ) 



d 



a 3 dt 



Eqs. d2H, fl2J2QD, and (JET) give, 



(3.1) 



(3.2) 



5p-R 5p - - R 2 5X - R 3 5(f> = 

k 2 f 2A; 2 



i?o = -3i7(l + c: 



6# a 2 /t 2 



- 3{l + w)p 



i?2 = 



2A: 4 , / /c 2 

3(1 + 4 



3a 4 #K 2 
/t 2 f 2A; 2 



2H > 



3a 2 H 



6H 1 K 2 a 2 



3(l + u>)Pf ( K 



XX 



k 2 k 2 (9k 2 

R 3 = --<PK X + _ { -£j - 3(1 + ti,)p ^ ~ <P 2 K 



(3.3) 
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where we have used 8p = c 2 5p. Therefore, from Eqs. (13.ip — (13.3p . we find 



d 
dt 



( Sp\ ( Ro Ri R 2 R3 \ 

* P Pi P2 Ps 

SX Eq E± E2 E3 

\ 5<j> J \ -0 1/0 ) 



5X 



We rewrite the above equations to the dimensionless forms for the later convenience: 





( Sp/p \ 


( A n 


A 12 


Ai 3 


A u \ 




d 




A 2 i 


A22 


A23 


A 2 4 




dN 


6X/X 


A 31 


A 32 


A33 


A34 


5X/X 




\ s<t>/<f> J 


\ Au 


A 42 


A43 


A44 j 


V 50/0 / 



Here, Ay = 1,2,3,4) are given by 



An =3(w — c 2 



2k 2 



6H 2 a 2 K 2 



3(1 +w)p 



A12 

A13 

A u 
A 2 i 



2k 4 



3a 4 H 2 n 2 p 



+ 3{l + w) 1 + 



k 2 \ 
3a 2 H 2 ) 



2k 2 



12H 2 \n 2 a 2 p 
k 2 1 



3(1 + w] 



K 



x 



xx 



2 -<Jxj>(K, x + j> 2 K xx ) - (1 + w)—^^ - 2 i^) 

K 2 p 

OH 2 



A 



22 



3a 2 H 2 



A 23 = - 
A31- 

A32 



12H 2 



K. x + <b 2 K, 



XX 



A24 = 



K 



6H 2 



K,-6 2 K 



,x<t> ' 



K 2 pK t) 



A 33 = -2 



H 2 {K. X + 4> 2 K. XX ) 
2k 2 K, x 20 
''a 2 H 2 K )X + <p 2 K )XX ~ H~4> 

d_ 
dN 



+ 



K 



#0 2H 2 



{K tX + <p 2 K tXX ] 



A34 = ~ 2^ 



K 



<P(K X + <P 2 K XX )\2H 2 



+ 



1 d 



a 3 dN 
An =0, A A2 = 



a 4>K jX <p 



H<f) 



A 43 = 



2H<p 



A 



li 
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B. Evolution of the matter perturbation 



In the last subsection, we have derived a simultaneous differential equations for the per- 
turbative quantities. We will decompose the equations to give a single differential equation 
of the matter perturbation. As a first step, we consider the case without the scalar field. 
Then, we obtain following equation from the first and the second lines of Eq. ( 13. 5p . 



d 2 S / A d , . „ \ dS 



A n A 22 - A 12 A 21 - ^ + A n -^ In |A 12 |^ 5 . 



(3.6) 



This is the differential equation for the matter density perturbation without approximation 



in the ACDM model [20|]. If we use the sound speed c s considered subhorizon approximation, 
that is, not d 6 ~ H5 but H 2 6 <C k 2 6/a 2 , Eq. (13. 6 j) is written as, 

d 2 5 fl a 3 ^ffk 2 d5 

+ {^h + 3{w ~ Cs2)(1 + 3w+ 3Wcs) ~ l {1 + 3w)(1 + Wcff) + (b^)" 1 ) } 5 • 

(3.7) 

Here, we have used w e s = —2H/3H 2 — 1 and (1 + w)k 2 p/ (3H 2 ) = 1 + w e g. Now we consider 
the case with the fc-essence field <fi. From Eq. (13. 5ft . we obtain the fourth order differential 
equation, 

d A 5 ,d 3 S ,d 2 5 d5 



where the definition of the coefficients M 3 , M 2 , Mi, and M are written in Appendix A. 
The expressions of the coefficients in the differential equation (13. 8p are so complicated that 
it is difficult to analyze the differential equation without any approximations. If we use 
small scale approximation, k 2 /a 2 3> H 2 , we can only consider the leading order terms in the 
expansion by the power series of a 2 H 2 /k 2 in the coefficients in Eq. (13.81) as follows: 



M. 



7 ^o 2 R 20 k 2 P , 3k 2 2 A> 6i*T x 

7 + 3c B — 6w r — 3(1 + w) — — 1 



H<j> H 2 H 2 4K >X + (j) 2 K,xx \ Hcj) 

<j>(10K x + 2 A> x )(0A> x + K t4>x ) 30A> 



HK tX (4K tX + 2 A> X ) H4>(K X + 2 A> X )(4A> + 2 A>x) 

x {3(1 + w)k 2 P K, x - 3k 2 <P 2 K 2 x + 2K M - 2<j)K^ x - 6H<j)K^ x - 2$K MX 




a 2 H 2 
+ 



0(1) 



(1 + w)k 2 P K 2 <p 2 K yX 



(3.9) 
(3.10) 



2H 2 



2H 2 




K 



x 



Qw 



(1 + w)k 2 P K 2 <j) 2 K iX 



K !X + (j) 2 K jXX 
-2K. X + 6 2 K. XX - 64jK. 



+ ct 



H4> 



2H 2 2H 2 

-0(lOif x + <j> 2 K xx )(<pK xx + K t vc) 



+ ci 



AK X + ct> 2 K. xx 
3K X 



- c a 



HK, X (4K, X + <P 2 K, XX ] 



H<P(K, X + <P 2 K, XX )(AK. 



x 



9K 



XX) 



-3H<P 2 (4>K, XX + 2K 



<px) 



+ 3(1 + w)n 2 p(j)K )X - 3n 2 (p 3 K 2 x - 2(fxf)K >tj>x + 2(j>K. 



<t>xx 



<j}(j>X 



k 4 



c 2 K 



x 



a 4 H 4 \K x + (b 2 K 



O 



,xx 



0(1) 



a 2 H 2 



(3-11) 



(3.12) 



Here, we have deleted H, p, and the higher derivatives of them by using Eqs. ( 12. 5ft and ( 12.71) . 
Furthermore, we have used the field equation (12. 6|) to delete and also used derivatives 
of Eq. (12. 6 p to delete K xxx , K xxxx , and K xxxxx . Thus, these coefficients (I3.9p -( 13~T2|) 
can be expressed in different forms by using the background equations (l2.4p -( l2~7Tj) and their 
derivatives. 

If we use the WKB approximation under the condition, k 2 /a 2 3> H 2 , we can obtain the 
solution of Eq. (13.81) as follows, 



S(N) 



4 

E 

i=l 



Ci(N)exp 



N 



an 



(3.13) 



Here, Cj(iV)'s are functions determined by the leading terms in Eqs. (13 . 9 1) ( 13TTT]) as the terms 
proportional to k 3 /(aH) 3 in Eq. (13.81) . and Aj's are determined by the terms proportional 
to (k/aH) 4 as Ai = — ic s , A 2 = ic s , A 3 = — ic^, and A 4 = ic^, respectively. The sound speed 
in fc-essence model = {pcj>),x/ {Pcj>),x S is given by K tX /(K^ x + <p 2 K >xx ), where p^ and 
Ptf, are the energy density and pressure of the fc-essence scalar field, respectively. If c s is a 
real number and c^((p, X) is a real function, then the solution ( 13 . 1 3[) oscillates rapidly in the 
small scale k/a^> H since Aj's are pure imaginary. We will investigate this oscillation later. 
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If we require the solution (I3.13P to be a real function, we should replace it with the following 
form: 



N 



S(N) =d(N) cos 

+ C 3 (N) cos 



an 



N 



N dN'^ 
an 



+ C 2 (N) sin 
+ C 4 (N) sin 



dN'^- 
an 



N 



dN 



i c<f>k 
~a~H 



(3.14) 



where d(N), C 2 (N), C 3 (N), and C 4 (N) are real functions of N. d(N) and d{N) are 
expressed as follows by to solve the third order equation of k/(aH) of Eq. (13.81) . 



jfih\Ci(N)\ = ^ln\d(N)\ = ~(-l + 6™ 



3c 



(3.15) 



The above expression gives dln\Ci(N)\/dN = din \d(N)\/dN < 0, therefore, these are 
not growing modes of the matter density perturbation. The behavior of these solutions 



Ci(iV)cos f N dN'c s k/(aH) and C 2 (iV) sin f lv dN f c s k/(aH) 
havior of the solutions of Eq. (13.71) when c s ^ 0. In the following, we will focus on the matter 
dominant era and the subsequent era, then the equation of state for matter and the sound 
speed of matter can be treated as w = and c s = 0, respectively. The term proportional 
to /c 4 /(a 4 if 4 ) in Eq. ( 13. 121) . then, vanishes, and we should evaluate more lower order terms. 
When c 2 = w = 0, M is expressed as 

k 2 3 



-N 



are completely the same be- 



M 0" 4^2 2 



n m 5 + o(i) 



(3.16) 



For Mq in fl3~T6|) . Ci(iV)cos J N dN'c s k/{aH) and C 2 (iV) sin J N dN'c s k/(aH) 



the solutions of Eq. 



are not 

In the following subsections, we investigate the solutions and 



estimate the time evolution of C^(N) and C^N) in case of quintessence model and general 
fc-essence model. 

C. Quintessence model 

When we restrict the fc-essence action to be the action of quintessence model (K((j), X) = 
—X + V((f))), the forms of functions Cs(N) and C^{N) are given by 



(1 , ,_„ rtl d , 3 d> 6 6 

— In \C 3 (N)\ = — ]n\CAN)\ = — + — -r 

dN 1 31 M dN 1 41 n 2 h<}> 2H4> 2#0 



3H + V, 



(3.17) 
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We can simplify this expression (13 . 1 7[) by using the equation of motion for the quintessence 
model, (p + 3H(p + = as follows: 

^\p, m = ^m\pm\ = ^. (3.18) 

Then we find Cs(N), C^N) oc <fi(N). The effective evolution of the growth factor 
dln5(N)/dN is determined by the evolution of f cS = din \C 3 (N)\/dN = din \C 4 (N)\/dN = 
4>/(H(j)) because the exponential factor of the solution ( 13. 14ft is oscillating. If we consider 
a quintessence model with an exponential potential, which is V(4>) = B 2 e~^/^ 2K ^ , we know 



the exact solutions 21] 



a 3 (t) =(uit + u 2 ) \^uxK 2 B 2 t 3 + ^u 2 K 2 B 2 t 2 + v x t + v^j , (3.19) 

, m ___y^2_ ] uit + u 2 . . 

m ~ V3k \u x k 2 B 2 V + \u 2 k 2 BH 2 + v x t + v 2 ' [6 - ZU) 

where u x , u 2 , Vi, and v 2 are arbitrary constants of integration. We now impose u 2 = v 2 = 

to satisfy the conditions a(0) = and f^(0) = because we do not consider the inflationary 

regime here. Then the effective growth factor f e g = <p/(H(f)) is given by 

'-M—t + L^ + a* - (3 ' 21) 

Here, the constants u± and V\ are restricted to satisfy po = 4u 1 ^ 1 /3k 2 by Eq. (12.4)) when the 
present scale factor a is chosen to be unity, so that u x v x > because p > 0. Therefore, 
fes(t) is a monotonically decreasing function and the value of f e g is restricted to be —3/4 < 
fcS < 3/2 because lim^o /efi = 3/2 and lim^oo / e fr = —3/4. The current value of f e s is 
estimated by using the parameter t s = 15.8 x 10 9 years defined by t 2 = Avi/k 2 B 2 ui. The 



value of t s follows from the conditions f2 m o = 0.3 and H = 70kms 1 Mpc 1 21]. With the 



present time to = 13 x 10 9 years, we find f c s(t ) « 0.21, which does not contradict with the 



observational results 



22] 



On the other hand, the remaining two solutions in four linearly independent solutions are 
determined as follows. In quintessence model, the coefficient M' Q is expressed as 

< = - ^1*1*8 + 0{1) • (3.22) 

By taking account of the second order terms of k/(aH) of coefficients of Eq. (j3.8|) . we find 
the following equation, 

d 2 5 (\ 3 \ d5 3^ P , 

o " O^ff ^ ~ = 0, (3.23) 



dN 2 V 2 2 "7 dN 2 
12 



which will yield quasi-static solutions. This equation ( I3.23|) is same as Eq. (I2.24|) . In case 
of exponential potential quintessence model we have considered, solutions of Eq. f!3.23j) are 
sub-leading ones because the growth rate of them are less than 3/2 in matter dominant era. 



D. General fc-essence model 



In case of general /c-essence model, Eq. (13.231) is written as: 



f d 2 5 (\ 3 \ d6 3 



This equation fl3T24"j) equals to Eq. (]2]Z1D when c ^ i.e. K >x ^ 0. 

On the other hand, the coefficient of the oscillating solution C (N) is determined as 

_ 5 QK+x + #X,xx 0(10^ x + <j> 2 K >xx )(]>K xx + K,ix) 



H<p 4 HK tX 2HK tX (4,K tX + (j) 2 K >xx 

5 3K. X 



4H(K tX + (^K )XX ) E{K )X + (P 2 K >xx )(4K tX + 4> 2 K XXi 
x (j>K i4JC + 3^K XX + 3 /^, xx + 3 <^,xxx) • (3.25) 

Here, the derivative of the field equation (I2.6P is also used to simplify the equation. 
We will discuss the possibility of "reconstruction" of the effective growth factor / cff = 
d\n\Cs(N)\/dN = c/ln |C4(iV)|/<i./V in the following. Before considering the reconstruction 
of the effective growth factor f C R , we review on the reconstruction of background geometry. 
The reconstruction, we call here, is the method for constructing the action which yields 
some given behavior of quantities, e.g., a(t), H(t), or f(N) = d\n5(N)/dN. We can always 
express <p(t) as <p(t) = g(t) by using some bijection g(t) because the redefinition of <fi{t) can 
be absorbed into the redefinition of the function K((p,X). If we choose g(t) as = t, the 
action of /c-essence model is reconstructed for arbitrary growth of the background space-time 



(for an arbitrary a(t) or H(t)) as follows [23j : 



{</>) = [wp+ — + — 













V 


3H 2 


2H 


K 2 
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H 2 



(3.26) 



The terms which are proportional to the square or higher power of (X — 1/2) do not con- 
tribute to the background evolution of space-time. Therefore, we can choose their coefficients 
K^((f>) with (n > 1) as arbitrary functions of 0. 

We now consider that the possibility of reconstruction of the effective growth factor by 
using the arbitrariness of K^((p). First, we should note that Eq. (I3.25P is rewritten as 
follows: 



1 d , 
H In 

AdN 



+ - dN 



4-ki\4Kr + j, 2 K :XX \, (3.27) 



where we have used C(N) oc Cs(N) oc C±(N) as shown in (13.251) . From the above equation 
(I3.27p . we immediately obtain the following relation: 

C{N) 4 K% = const, x 4 (X X + j?K >xx ){AK >x + <p 2 K, xx ) 4 . (3.28) 

The reconstruction of the effective growth factor f eS is completed if f e s is integrable and 
Eq. (I3.28P is algebraically solvable for K xx because C(N) is expressed by using the in- 
tegration of f c f[. However, we cannot complete the reconstruction because Eq. (I3.28|) is a 
quintic algebraic equation with respect to K xx and hence it is not generally solvable. It 
should be noted that it is possible to execute the reconstruction approximately. By redefin- 
ing a{N) = <p 2 K, xx /K <x and y(N) = C 4 (X)/(const. x <p 4 K 2 x ), we can rewrite Eq. (13T28D 
as 

y(N) = (l + a(iV))(4 + «(iV)) 4 . (3.29) 

Thus, we find that there are one or more real solutions for a. To be concrete, in case 
y < — 2 8 3 5 /5 5 or < y there is one real solution for Eq. (13.291) . in case y = — 2 8 3 5 /5 5 ,0 
there are two real solutions, and in case — 2 8 3 5 /5 5 < y < there are three real solutions. So 
that it is possible to define a(N) for the arbitrary behavior of y(N). It should be, however, 
noted that the function a(N) must be a multi- valued function or a discontinuous function 
if we consider the continuous function y(N) which crosses over the region [— 2 8 3 5 /5 5 , 0]. 
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E. Oscillation 



Since we find the ratio between Cs(N) and C^{N) is a constant factor from Eq. (13.251) . 
we can rewrite the oscillating solutions in the following form: 



where u is an arbitrary real constant. The frequency of the solution ( 13.301) is given by 
ck/(a H ) » 300 when we use the values, H = h x (9.777752 Gyr)" 1 and k = 0.1/iMpc -1 . 
This corresponds to the frequency around the present time in quintessence model because 
C0 = c is satisfied in quintessence model. We now examine how often the matter density 
perturbation becomes 0. We set 300Az « ir and obtain Az ~ 7r/300 because function sin^ 
vanishes in the period of 9 by tt. There are points that the matter density perturbation 
vanishes, for the period 0.01 of the redshift z. Generally, the frequency of the solution (I3.30f) 
is given by ck(l + z')/(aoH(z = z')) when we consider it around z = z'. The difference 
of N, AN(z') = N(z' - Az) - N(z'), is expressed as ln(l + z')/(l + z' - Az) due to the 
definition of N. Hence, AN « Az/(1 + z') is justified when |1 + z'\ 3> |Az|. The frequency 
corresponding to z, then, is expressed as ckjia.^Hi^z = z')). The difference of the frequency 
corresponding to z between around z = and around z = z' is only H(z). If we use a 
scale factor of the ACDM model, a(t) = (p m o/ Pa) 1 ^ 3 sinh 2//3 [-\/3fi; 2 pA/4 t], we can evaluate 
the concrete value of the frequency because the ratio H(z = z')/H(z = 0) is determined. 
In case of z = 1 and k = 0.01/i Mpc -1 , the frequency is approximately one seventeenth 
times of that in case of z = and k = 0.1/iMpc" 1 because the ratio H(z = 1)/H(z = 0) is 
approximately determined 1.7. Here, we have used Q\ = 0.74, Q m o = 0.26, and h = 0.72. 
Thus, the distance of the points that the matter density perturbation vanishes becomes 
bigger than that in case of z = and k = 0.1/zMpc" 1 , and the points appears in the interval 
of 0.18 in the redshift. 

We now consider the variation of the wave number k with a fixed redshift. The fre- 
quency corresponding to k, then, is expressed by J N dN'c ( / } /(aH). It can be replaced with 
J* dtc^ I 'a(t) if we use dN = H(t)dt. Here, we assume the bottom point of the domain 
of the integration as a time of the radiation-matter equality t cq . We also investigate it 
in case of = c and a(t) = (p m o/ Pa) 1 ^ 3 sinh 2 / 3 [A/3ft 2 pA/4 t] in order to evaluate the 
value of the frequency corresponding to k. The present value of the frequency is given 




(3.30) 
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by c J^dt/a(t) » 1.45 x 10 4 Mpc, where we have used t„ = 13.7 x 10 9 years, t eq = 4000 
years, and the sidereal year equals 31558149.8 s. So that, the amount of the variation of k 
in order to make a phase shift with ir is A A; » 2.17 x 10~ 4 Mpc _1 . It is so small that the 
observational verification should be very difficult. On the other hand, A/c in z = 1, which 
is given by setting the upper end of the domain of the integration to be t = 6.0 x 10 9 years, 
is approximately 1.3 times bigger than that in z = 0, but it is also very small. 

IV. CONCLUSIONS AND DISCUSSIONS 

We have estimated the time evolution of the matter density perturbation without using 
subhorizon approximation but using only small scale approximation k/(aH) 3> 1 in this 
paper. As a result, the characteristic behaviors of the matter density perturbation, the 
oscillation caused by the sound speed in /c-essence model, and the effective growth factor 
depends on K t x, K y xx, 4>, and their derivatives, have been discovered. From a viewpoint 
of verification of dark energy, the effective growth factor is not very important because 
even the simple example of quintessence model does not induce a great difference from the 
ACDM model, and the possibility of approximate reconstruction of the effective growth 
factor has been shown. However, the oscillation caused by the sound speed in /c-essence 
model is a peculiar behavior which is not explained by the ACDM model, so that we can 
find a clear difference in it. As seen in the last section, the oscillation could be, moreover, 
observable if we fix the wavenumber k and move the redshift z around k = 0.01/iMpc -1 
and z — 1. On the other hand, second order differentiated equation of the matter density 
perturbation which will give quasi-static solutions has been found, too. This differential 
equation is same as that, for which the subhorizon approximation have been applied when 
Ctj, 7^ 0. Therefore, there is a possibility that the quasi-static solutions have a growing mode 
but the oscillating solutions do not have a growing mode, and thus the oscillation could not 
be able to be observed. Several conditions, however, will be imposed on the action to realize 
such a situation. Consequently, if such an oscillation is not observed, /c-essence model would 
be severely constrained. On the other hand, if such an oscillation is observed, the existence 
of the scalar field dark energy such as /c-essence model or the dark matter given by a scalar 
field may be verified. 
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Appendix A: Derivation of the forth order differential equation 

The fourth order differential equation (13.81) . 

A d 3 5 A d 2 6 n d5 „ r r n 
— + M 3 — + M 2 — + ^- + ^£ = 0, (Al) 

which is given by the Gaussian elimination like procedure: 

1. Differentiating the first line of Eq. ( 13. 5ft with respect to N. 

2. Eliminating the terms proportional to dfy/dN, d(5X/X)/dN, and d(5<fr/(j))/dN from 
the equation derived in the first step by using the second, third, and fourth lines of 
Eq. (1331). 

3. Differentiating the equation given in the second step with respect to N. 

4. Eliminating the terms proportional to dfy/dN, d(5X/X)/dN, and d(5<fr/(j))/dN from 
the equation derived in the third step by using the second, third, and fourth lines of 
Eq. (JS3J. 

5. Differentiating the equation given in the fourth step with respect to N. 

6. Eliminating the terms proportional to d^/dN, d(5X/X)/dN, and d(5<fr/ '<fi) /dN from 
the equation derived in the fifth step by using the second, third, and fourth lines of 
Eq. (13. 5p . Then we obtain the fourth order differential equation for 5 including the 
extra quantities SX/X, and 



7. By using the equations obtained in the second and fourth steps and the first line of 
Eq. (I3.5p . we can delete the terms proportional to 5X/X, and 5<j)/<j) in the equation 
obtained in the sixth step, and then we find Eq. (1A1|) . 
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Eq. ( 1Aip is the differential equation for the matter density perturbation in fc-essence 
without any approximation. The coefficients in flAlj) are defined by 

4 

M 3 



gL4i 7 



Mo 



-A n + D(a t P 7 i)AiJA lf) 

4 dA 4 \ 

+ AuA^ — Ai~/—j-j^. — Ai y AiiAip J , 

i=2 i=2 ' 

.4 4 4 . 



dN 



d 2 A^ 



dN 2 



a=2 Q ,^,7=2&a^/3^7 

+ a 11 a 1 ^^ + A ^ A w E A ^ + E ^w Ad ^ 

d=2 d=2 
4 ^ 4 

+ A ^^- + A W E A ^ A edA dl 

d=2 d,e=2 

d?Aip dA\p s—^ 

~ 7 dN 2 ~ ll ~dN~ ~ 7 ^ 

d=2 

4 ^ 4 ^ 4 x 

- 2A ^ E - E - A ^ E ^^a^ ) 

d=2 d=2 d,e=2 ' 

, , d 2 A-[-[ v — -\ rfyliry , v — -\ , dA n -\ v — -\ 

UlS - 3 *v^ - 3 E w- 4 - - 2 E ^« -w ~ E 

a=2 a=2 a„S=2 



D(a, p,<y)H 2A W -^ + Ay J2 A ^ A di 

a,/3,7=2&a^/3^7 ^ ^ d=2 



dA^ 

] lLN~ 



+ 2 E + E A ^ d ~dw + E ^-^) ( A 

d=2 d=2 d,e=2 ' ^ 

4 4 x 

+ Ay A i d A dl - A ii—^f- ~ A n E A uAdp J 

d=2 d=2 ' 

v d=l 7 v d=2 

+ Ay A ™~dN + ^ E A ^ A ed A d, - A l,^f 
d=2 d,e=2 

- 2A ^ E ^dW Ad ? - A ^< E A ^^ - Ai ^ E A ^ A edA dp \ } , 

d=2 d=2 d,e=2 ' ' 

_ d?A\\ sr-^ d 2 Ai a , dAi a dA a i ^-^ d 2 A a i 



dN 3 ^ dN 2 U1 ^ dN dN ^ 1U dN 2 

a=2 a=2 a=2 
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3 E ^j^ A(3aAai ~ 2 E Aipd ^f Aai ~ E Ai ^ a ^dW 

a,/3=2 a,P=2 a,/3=2 



a., 

- A ly A^Ap a A al + jD ^' /3,7 ^{( Al1 ^^ 

+ E d ~w A v + A » E aJ^- + Mi E M,M d A d , 

d=2 d=2 d,e=2 



dA dl 
ld dN 



d=2 d=2 

- A w A ieA ed A dl ) Uip^ + A 1P A ^ A h ~ A ^^dW 

d,e=2 ' ^ f=2 

- A„ ± A lfAfe ) - U,^? 4- 2^ £ ^* 

f=2 / V d=2 



" A V E + ^ E ^4*4* " ^7^T 

d=2 d,e=2 

4,. 4 J /I 4 > 

^ — \ u^x^ ^ — \ dA d p ^ — \ 

d=2 d=2 d,e=2 y 



x A n 



^ + Au E - V^F - ^ E } , (A5) 

f=2 f=2 J ) 

4 , A 4 , . 



k v d=2 d=2 

\ — V 4 4 4 4 (P A\ a ^ — \ dAld 

+ A la A le A ed A dp - A xf} - 2A 1/3 ~^Jf Ada 

d,e=2 d=2 

~ ^ E M^ - M„ t, M e A cd A da ) U/-^L 

d=2 d,e=2 ' ^ 

+ a w A v A h - A ^^r ~ Mi E A v A fp) 

f=2 f=2 J 

-(aJ& + 2A V ± d 4i Ail + A„ E Aj A - 



d=2 d=2 



dN 



a \ — v 444 a A\fi ^ — \ dA\ d 

+ A W 2^ A le A ed A dl - Airy—jj^- - 2Ai 7 2^ —[^ A df} 

d,e=2 d=2 

- ^X>«^r - All ± A le A cd A dB ) UJ^. 

d=2 d,e=2 ' ^ 
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dAi a 



+ A la J2 AifAfp ~ A ^~^ ~ A W Yl Al f A 



x 



/=2 

d 3 A la 
dN 3 
1 



fa 



/=2 



\^ d 2 Ai d ^ dA^dAda . d 2 A da 

6 2^^w Ad « + 6 2^ dN dN + 2^ A ™-dw 

d=2 d=2 d=2 



dA f 



dA 



^ dN ^ dN ^ dN 



d,e=2 
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